Genomics Research of Lifetime Depression in the Netherlands: The BIObanks Netherlands Internet Collaboration (BIONIC) Project.

In this cohort profile article we describe the lifetime major depressive disorder (MDD) database that has been established as part of the BIObanks Netherlands Internet Collaboration (BIONIC). Across the Netherlands we collected data on Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition (DSM-5) lifetime MDD diagnosis in 132,850 Dutch individuals. Currently, N = 66,684 of these also have genomewide single nucleotide polymorphism (SNP) data. We initiated this project because the complex genetic basis of MDD requires large population-wide studies with uniform in-depth phenotyping. For standardized phenotyping we developed the LIDAS (LIfetime Depression Assessment Survey), which then was used to measure MDD in 11 Dutch cohorts. Data from these cohorts were combined with diagnostic interview depression data from 5 clinical cohorts to create a dataset of N = 29,650 lifetime MDD cases (22%) meeting DSM-5 criteria and 94,300 screened controls. In addition, genomewide genotype data from the cohorts were assembled into a genomewide association study (GWAS) dataset of N = 66,684 Dutch individuals (25.3% cases). Phenotype data include DSM-5-based MDD diagnoses, sociodemographic variables, information on lifestyle and BMI, characteristics of depressive symptoms and episodes, and psychiatric diagnosis and treatment history. We describe the establishment and harmonization of the BIONIC phenotype and GWAS datasets and provide an overview of the available information and sample characteristics. Our next step is the GWAS of lifetime MDD in the Netherlands, with future plans including fine-grained genetic analyses of depression characteristics, international collaborations and multi-omics studies.

Genetic confounding in the association of early motor development with childhood and adolescent exercise behavior.

INTRODUCTION: Early motor development has been found to be a predictor of exercise behavior in children and adolescents, but whether this reflects a causal effect or confounding by genetic or shared environmental factors remains to be established. METHODS: For 20,911 complete twin pairs from the Netherlands Twin Register a motor development score was obtained from maternal reports on the timing of five motor milestones. During a 12-year follow-up, subsamples of the mothers reported on the twins' ability to perform seven gross motor skills ability (N = 17,189 pairs), and weekly minutes of total metabolic equivalents of task (MET) spent on sports and exercise activities at age 7 (N = 3632 pairs), age 10 (N = 3735 pairs), age 12 (N = 7043 pairs), and age 14 (N = 3990 pairs). Multivariate phenotypic and genetic regression analyses were used to establish the predictive strength of the two motor development traits for future exercise behavior, the contribution of genetic and shared environmental factors to the variance in all traits, and the contribution of familial confounding to the phenotypic prediction. RESULTS: Significant heritability (h2) and shared environmental (c2) effects were found for early motor development in boys and girls (h2 = 43-65%; c2 = 16-48%). For exercise behavior, genetic influences increased with age (boys: h2age7 = 22% to h2age14 = 51%; girls: h2age7 = 3% to h2age14 = 18%) paired to a parallel decrease in the influence of the shared environment (boys: c2age7 = 68% to c2age14 = 19%; girls: c2age7 = 80% to c2age14 = 48%). Early motor development explained 4.3% (p < 0.001) of the variance in future exercise behavior in boys but only 1.9% (p < 0.001) in girls. If the effect in boys was due to a causal effect of motor development on exercise behavior, all of the factors influencing motor development would, through the causal chain, also influence future exercise behavior. Instead, only the genetic parts of the regression of exercise behavior on motor development were significant. Shared and unique environmental parts of the regression were largely non-significant, which is at odds with the causal hypothesis. CONCLUSION: No support was found for a direct causal effect in the association between rapid early motor development on future exercise behavior. In boys, early motor development appears to be an expression of the same genetic factors that underlie the heritability of childhood and early adolescent exercise behavior.

Using Instrumental Variables to Measure Causation over Time in Cross-Lagged Panel Models.

Cross-lagged panel models (CLPMs) are commonly used to estimate causal influences between two variables with repeated assessments. The lagged effects in a CLPM depend on the time interval between assessments, eventually becoming undetectable at longer intervals. To address this limitation, we incorporate instrumental variables (IVs) into the CLPM with two study waves and two variables. Doing so enables estimation of both the lagged (i.e., "distal") effects and the bidirectional cross-sectional (i.e., "proximal") effects at each wave. The distal effects reflect Granger-causal influences across time, which decay with increasing time intervals. The proximal effects capture causal influences that accrue over time and can help infer causality when the distal effects become undetectable at longer intervals. Significant proximal effects, with a negligible distal effect, would imply that the time interval is too long to estimate a lagged effect at that time interval using the standard CLPM. Through simulations and an empirical application, we demonstrate the impact of time intervals on causal inference in the CLPM and present modeling strategies to detect causal influences regardless of the time interval in a study. Furthermore, to motivate empirical applications of the proposed model, we highlight the utility and limitations of using genetic variables as IVs in large-scale panel studies.

Beyond the Global Brain Differences: Intraindividual Variability Differences in 1q21.1 Distal and 15q11.2 BP1-BP2 Deletion Carriers.

BACKGROUND: Carriers of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants exhibit regional and global brain differences compared with noncarriers. However, interpreting regional differences is challenging if a global difference drives the regional brain differences. Intraindividual variability measures can be used to test for regional differences beyond global differences in brain structure. METHODS: Magnetic resonance imaging data were used to obtain regional brain values for 1q21.1 distal deletion (n = 30) and duplication (n = 27) and 15q11.2 BP1-BP2 deletion (n = 170) and duplication (n = 243) carriers and matched noncarriers (n = 2350). Regional intra-deviation scores, i.e., the standardized difference between an individual's regional difference and global difference, were used to test for regional differences that diverge from the global difference. RESULTS: For the 1q21.1 distal deletion carriers, cortical surface area for regions in the medial visual cortex, posterior cingulate, and temporal pole differed less and regions in the prefrontal and superior temporal cortex differed more than the global difference in cortical surface area. For the 15q11.2 BP1-BP2 deletion carriers, cortical thickness in regions in the medial visual cortex, auditory cortex, and temporal pole differed less and the prefrontal and somatosensory cortex differed more than the global difference in cortical thickness. CONCLUSIONS: We find evidence for regional effects beyond differences in global brain measures in 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants. The results provide new insight into brain profiling of the 1q21.1 distal and 15q11.2 BP1-BP2 copy number variants, with the potential to increase understanding of the mechanisms involved in altered neurodevelopment.

Contributions of amyloid beta and cerebral small vessel disease in clinical decline.

INTRODUCTION: We assessed whether co-morbid small vessel disease (SVD) has clinical predictive value in preclinical or prodromal Alzheimer's disease. METHODS: In 1090 non-demented participants (65.4 ± 10.7 years) SVD was assessed with magnetic resonance imaging and amyloid beta (Aß) with lumbar puncture and/or positron emission tomography scan (mean follow-up for cognitive function 3.1 ± 2.4 years). RESULTS: Thirty-nine percent had neither Aß nor SVD (A-V-), 21% had SVD only (A-V+), 23% Aß only (A+V-), and 17% had both (A+V+). Pooled cohort linear mixed model analyses demonstrated that compared to A-V- (reference), A+V- had a faster rate of cognitive decline. Co-morbid SVD (A+V+) did not further increase rate of decline. Cox regression showed that dementia risk was modestly increased in A-V+ (hazard ratio [95% confidence interval: 1.8 [1.0-3.2]) and most strongly in A+ groups. Also, mortality risk was increased in A+ groups. DISCUSSION: In non-demented persons Aß was predictive of cognitive decline, dementia, and mortality. SVD modestly predicts dementia in A-, but did not increase deleterious effects in A+. HIGHLIGHTS: Amyloid beta (Aß; A) was predictive for cognitive decline, dementia, and mortality. Small vessel disease (SVD) had no additional deleterious effects in A+. SVD modestly predicted dementia in A-. Aß should be assessed even when magnetic resonance imaging indicates vascular cognitive impairment.

Longitudinal multi-omics study reveals common etiology underlying association between plasma proteome and BMI trajectories in adolescent and young adult twins.

BACKGROUND: The influence of genetics and environment on the association of the plasma proteome with body mass index (BMI) and changes in BMI remains underexplored, and the links to other omics in these associations remain to be investigated. We characterized protein-BMI trajectory associations in adolescents and adults and how these connect to other omics layers. METHODS: Our study included two cohorts of longitudinally followed twins: FinnTwin12 (N = 651) and the Netherlands Twin Register (NTR) (N = 665). Follow-up comprised 4 BMI measurements over approximately 6 (NTR: 23-27 years old) to 10 years (FinnTwin12: 12-22 years old), with omics data collected at the last BMI measurement. BMI changes were calculated in latent growth curve models. Mixed-effects models were used to quantify the associations between the abundance of 439 plasma proteins with BMI at blood sampling and changes in BMI. In FinnTwin12, the sources of genetic and environmental variation underlying the protein abundances were quantified by twin models, as were the associations of proteins with BMI and BMI changes. In NTR, we investigated the association of gene expression of genes encoding proteins identified in FinnTwin12 with BMI and changes in BMI. We linked identified proteins and their coding genes to plasma metabolites and polygenic risk scores (PRS) applying mixed-effects models and correlation networks. RESULTS: We identified 66 and 14 proteins associated with BMI at blood sampling and changes in BMI, respectively. The average heritability of these proteins was 35%. Of the 66 BMI-protein associations, 43 and 12 showed genetic and environmental correlations, respectively, including 8 proteins showing both. Similarly, we observed 7 and 3 genetic and environmental correlations between changes in BMI and protein abundance, respectively. S100A8 gene expression was associated with BMI at blood sampling, and the PRG4 and CFI genes were associated with BMI changes. Proteins showed strong connections with metabolites and PRSs, but we observed no multi-omics connections among gene expression and other omics layers. CONCLUSIONS: Associations between the proteome and BMI trajectories are characterized by shared genetic, environmental, and metabolic etiologies. We observed few gene-protein pairs associated with BMI or changes in BMI at the proteome and transcriptome levels.

Phenotypic but not genetically predicted heart rate variability associated with all-cause mortality.

Low heart rate variability (HRV) has been widely reported as a predictor for increased mortality. However, the molecular mechanisms are poorly understood. Therefore, this study aimed to identify novel genetic loci associated with HRV and assess the association of phenotypic HRV and genetically predicted HRV with mortality. In a GWAS of 46,075 European ancestry individuals from UK biobank, we identified 17 independent genome-wide significant genetic variants in 16 loci associated with HRV traits. Notably, eight of these loci (RNF220, GNB4, LINCR-002, KLHL3/HNRNPA0, CHRM2, KCNJ5, MED13L, and C160rf72) have not been reported previously. In a prospective phenotypic relationship between HRV and mortality during a median follow-up of seven years, individuals with lower HRV had higher risk of dying from any cause. Genetically predicted HRV, as determined by the genetic risk scores, was not associated with mortality. To the best of our knowledge, the findings provide novel biological insights into the mechanisms underlying HRV. These results also underline the role of the cardiac autonomic nervous system, as indexed by HRV, in predicting mortality.

Genome-Wide DNA Methylation Profiles in Whole-Blood and Buccal Samples-Cross-Sectional, Longitudinal, and across Platforms.

The field of DNA methylation research is rapidly evolving, focusing on disease and phenotype changes over time using methylation measurements from diverse tissue sources and multiple array platforms. Consequently, identifying the extent of longitudinal, inter-tissue, and inter-platform variation in DNA methylation is crucial for future advancement. DNA methylation was measured in 375 individuals, with 197 of those having 2 blood sample measurements ~10 years apart. Whole-blood samples were measured on Illumina Infinium 450K and EPIC methylation arrays, and buccal samples from a subset of 58 participants were measured on EPIC array. The data were analyzed with the aims to examine the correlation between methylation levels in longitudinal blood samples in 197 individuals, examine the correlation between methylation levels in the blood and buccal samples in 58 individuals, and examine the correlation between blood methylation profiles assessed on the EPIC and 450K arrays in 83 individuals. We identified 136,833, 7674, and 96,891 CpGs significantly and strongly correlated (>0.50) longitudinally, across blood and buccal samples as well as array platforms, respectively. A total of 3674 of these CpGs were shared across all three sets. Analysis of these shared CpGs identified previously found associations with aging, ancestry, and 7016 mQTLs as well.

Effects of smoking on genome-wide DNA methylation profiles: A study of discordant and concordant monozygotic twin pairs.

Background: Smoking-associated DNA methylation levels identified through epigenome-wide association studies (EWASs) are generally ascribed to smoking-reactive mechanisms, but the contribution of a shared genetic predisposition to smoking and DNA methylation levels is typically not accounted for. Methods: We exploited a strong within-family design, that is, the discordant monozygotic twin design, to study reactiveness of DNA methylation in blood cells to smoking and reversibility of methylation patterns upon quitting smoking. Illumina HumanMethylation450 BeadChip data were available for 769 monozygotic twin pairs (mean age = 36 years, range = 18-78, 70% female), including pairs discordant or concordant for current or former smoking. Results: In pairs discordant for current smoking, 13 differentially methylated CpGs were found between current smoking twins and their genetically identical co-twin who never smoked. Top sites include multiple CpGs in CACNA1D and GNG12, which encode subunits of a calcium voltage-gated channel and G protein, respectively. These proteins interact with the nicotinic acetylcholine receptor, suggesting that methylation levels at these CpGs might be reactive to nicotine exposure. All 13 CpGs have been previously associated with smoking in unrelated individuals and data from monozygotic pairs discordant for former smoking indicated that methylation patterns are to a large extent reversible upon smoking cessation. We further showed that differences in smoking level exposure for monozygotic twins who are both current smokers but differ in the number of cigarettes they smoke are reflected in their DNA methylation profiles. Conclusions: In conclusion, by analysing data from monozygotic twins, we robustly demonstrate that DNA methylation level in human blood cells is reactive to cigarette smoking. Funding: We acknowledge funding from the National Institute on Drug Abuse grant DA049867, the Netherlands Organization for Scientific Research (NWO): Biobanking and Biomolecular Research Infrastructure (BBMRI-NL, NWO 184.033.111) and the BBRMI-NL-financed BIOS Consortium (NWO 184.021.007), NWO Large Scale infrastructures X-Omics (184.034.019), Genotype/phenotype database for behaviour genetic and genetic epidemiological studies (ZonMw Middelgroot 911-09-032); Netherlands Twin Registry Repository: researching the interplay between genome and environment (NWO-Groot 480-15-001/674); the Avera Institute, Sioux Falls (USA), and the National Institutes of Health (NIH R01 HD042157-01A1, MH081802, Grand Opportunity grants 1RC2 MH089951 and 1RC2 MH089995); epigenetic data were generated at the Human Genomics Facility (HuGe-F) at ErasmusMC Rotterdam. Cotinine assaying was sponsored by the Neuroscience Campus Amsterdam. DIB acknowledges the Royal Netherlands Academy of Science Professor Award (PAH/6635).

The genetic information of people who smoke present distinctive characteristics. In particular, previous research has revealed differences in patterns of DNA methylation, a type of chemical modification that helps cells switch certain genes on or off. However, most of these studies could not establish for sure whether these changes were caused by smoking, predisposed individuals to smoke, or were driven by underlying genetic variation in the DNA sequence itself. To investigate this question, van Dongen et al. examined DNA methylation data from the blood cells of over 700 pairs of identical twins. These individuals share the exact same genetic information, making it possible to better evaluate the impact of lifestyle on DNA modifications. The analyses identified differences in methylation at 13 DNA locations in pairs of twins where one was a current smoker and their sibling had never smoked. Two of the genes code for proteins involved in the response to nicotine, the primary addictive chemical in cigarette smoke. The differences were smaller if one of the twins had stopped smoking, suggesting that quitting can help to reverse some of these changes. These findings confirm that DNA methylation in blood cells is influenced by cigarette smoke, which could help to better understand smoking-associated diseases. They also demonstrate how useful identical twins studies can be to identify methylation changes that are markers of lifestyle.

Longitudinal multi-omics study reveals common etiology underlying association between plasma proteome and BMI trajectories in adolescent and young adult twins.

Background: The influence of genetics and environment on the association of the plasma proteome with body mass index (BMI) and changes in BMI remain underexplored, and the links to other omics in these associations remain to be investigated. We characterized protein-BMI trajectory associations in adolescents and adults and how these connect to other omics layers. Methods: Our study included two cohorts of longitudinally followed twins: FinnTwin12 (N=651) and the Netherlands Twin Register (NTR) (N=665). Follow-up comprised four BMI measurements over approximately 6 (NTR: 23-27 years old) to 10 years (FinnTwin12: 12-22 years old), with omics data collected at the last BMI measurement. BMI changes were calculated using latent growth curve models. Mixed-effects models were used to quantify the associations between the abundance of 439 plasma proteins with BMI at blood sampling and changes in BMI. The sources of genetic and environmental variation underlying the protein abundances were quantified using twin models, as were the associations of proteins with BMI and BMI changes. In NTR, we investigated the association of gene expression of genes encoding proteins identified in FinnTwin12 with BMI and changes in BMI. We linked identified proteins and their coding genes to plasma metabolites and polygenic risk scores (PRS) using mixed-effect models and correlation networks. Results: We identified 66 and 14 proteins associated with BMI at blood sampling and changes in BMI, respectively. The average heritability of these proteins was 35%. Of the 66 BMI-protein associations, 43 and 12 showed genetic and environmental correlations, respectively, including 8 proteins showing both. Similarly, we observed 6 and 4 genetic and environmental correlations between changes in BMI and protein abundance, respectively. S100A8 gene expression was associated with BMI at blood sampling, and the PRG4 and CFI genes were associated with BMI changes. Proteins showed strong connections with many metabolites and PRSs, but we observed no multi-omics connections among gene expression and other omics layers. Conclusions: Associations between the proteome and BMI trajectories are characterized by shared genetic, environmental, and metabolic etiologies. We observed few gene-protein pairs associated with BMI or changes in BMI at the proteome and transcriptome levels.

Handedness and 23 Early Life Characteristics in 37,495 Dutch Twins.

In studies of singletons, a range of early-life characteristics have been reported to be associated with handedness, but some of these associations have failed to replicate. We examined associations between 23 early life characteristics with handedness in a large sample of 37,495 5-year-old twins. We considered three definitions of handedness: left-handedness (LH), mixed-handedness (MH), and non-right-handedness (NRH). Our main aim was to test whether the associations with sex, birth weight, gestational age, and season of birth - as reported in singletons - replicate in twins, and to examine twin-specific variables, including zygosity, chorionicity, birth order, and intertwin delivery time. Compared to previously published data from adults born as singletons (7.23%), the prevalence of NRH was higher in both twins (16.19%) and their parents (15.09%). In the twins, LH and NRH were associated with parents' LH. Male sex and lower gestational age were associated with NRH, and LH was associated with not being breastfed. MH was related to neurodevelopmental delays and higher externalizing problems later in childhood. Other previously reported associations were not replicated, and no twin-specific characteristics were related to handedness. These results emphasize the importance of considering multiple definitions of handedness and indicate a small number of replicated associations across studies.

Genetic and Environmental Effects on the Early Motor Development as a Function of Parental Educational Attainment.

INTRODUCTION: The contribution of genetic and environmental factors to individual differences in early motor development is still largely uncharted. This large-scale twin study establishes the genetic and environmental influences on the timing of motor milestones achievement, and it further tests whether the influences are moderated by parental education. METHODS: The twins came from families registered in the Netherlands Twin Register (NTR) from 1986 to 2016. In 30,256 complete twin pairs, mother-reported ages at which each twin was able to first-time roll from back to belly, sit unassisted, hands-and-knees crawl, stand up unaided, and walk independently were used to extract an early motor development factor. Parental education was dichotomized ("both parents with low/average education" vs "at least one parent with high education" with university degree as a threshold). RESULTS: Additive genetics explained 52% of the variance in motor development, the remaining 39% and 9% were explained by shared and nonshared environment separately. Mean age of achieving motor milestones tended to be higher in infants with high educated parents, and a moderation of parental education on the genetic and environmental variance in motor development was seen in female twins with larger heritability in the high educated parents group (64% vs 43%) paired to a lower shared environmental influence (28% vs 48%). Only 7%-8% of the variance was accounted for nonshared environmental factors, including measurement error. The pattern of results did not change when the degree of urbanicity, a correlate of parental education, was additionally considered. CONCLUSIONS: Genetic factors explain most of the individual differences in the timing of motor milestone achievement, but factors related to the shared home environment also play an important role in early motor development.

Genome-wide analysis identifies genetic effects on reproductive success and ongoing natural selection at the FADS locus.

Identifying genetic determinants of reproductive success may highlight mechanisms underlying fertility and identify alleles under present-day selection. Using data in 785,604 individuals of European ancestry, we identified 43 genomic loci associated with either number of children ever born (NEB) or childlessness. These loci span diverse aspects of reproductive biology, including puberty timing, age at first birth, sex hormone regulation, endometriosis and age at menopause. Missense variants in ARHGAP27 were associated with higher NEB but shorter reproductive lifespan, suggesting a trade-off at this locus between reproductive ageing and intensity. Other genes implicated by coding variants include PIK3IP1, ZFP82 and LRP4, and our results suggest a new role for the melanocortin 1 receptor (MC1R) in reproductive biology. As NEB is one component of evolutionary fitness, our identified associations indicate loci under present-day natural selection. Integration with data from historical selection scans highlighted an allele in the FADS1/2 gene locus that has been under selection for thousands of years and remains so today. Collectively, our findings demonstrate that a broad range of biological mechanisms contribute to reproductive success.

Genetically identical twin-pair difference models support the amyloid cascade hypothesis.

The amyloid cascade hypothesis has strongly impacted the Alzheimer's disease research agenda and clinical trial designs over the past decades, but precisely how amyloid-ß pathology initiates the aggregation of neocortical tau remains unclear. We cannot exclude the possibility of a shared upstream process driving both amyloid-ß and tau in an independent manner instead of there being a causal relationship between amyloid-ß and tau. Here, we tested the premise that if a causal relationship exists, then exposure should be associated with outcome both at the individual level as well as within identical twin-pairs, who are strongly matched on genetic, demographic and shared environmental background. Specifically, we tested associations between longitudinal amyloid-ß PET and cross-sectional tau PET, neurodegeneration and cognitive decline using genetically identical twin-pair difference models, which provide the unique opportunity of ruling out genetic and shared environmental effects as potential confounders in an association. We included 78 cognitively unimpaired identical twins with [18F]flutemetamol (amyloid-ß)-PET, [18F]flortaucipir (tau)-PET, MRI (hippocampal volume) and cognitive data (composite memory). Associations between each modality were tested at the individual level using generalized estimating equation models, and within identical twin-pairs using within-pair difference models. Mediation analyses were performed to test for directionality in the associations as suggested by the amyloid cascade hypothesis. At the individual level, we observed moderate-to-strong associations between amyloid-ß, tau, neurodegeneration and cognition. The within-pair difference models replicated results observed at the individual level with comparably strong effect sizes. Within-pair differences in amyloid-ß were strongly associated with within-pair differences in tau (ß = 0.68, P < 0.001), and moderately associated with within-pair differences in hippocampal volume (ß = -0.37, P = 0.03) and memory functioning (ß = -0.57, P < 0.001). Within-pair differences in tau were moderately associated with within-pair differences in hippocampal volume (ß = -0.53, P < 0.001) and strongly associated with within-pair differences in memory functioning (ß = -0.68, P < 0.001). Mediation analyses showed that of the total twin-difference effect of amyloid-ß on memory functioning, the proportion mediated through pathways including tau and hippocampal volume was 69.9%, which was largely attributable to the pathway leading from amyloid-ß to tau to memory functioning (proportion mediated, 51.6%). Our results indicate that associations between amyloid-ß, tau, neurodegeneration and cognition are unbiased by (genetic) confounding. Furthermore, effects of amyloid-ß on neurodegeneration and cognitive decline were fully mediated by tau. These novel findings in this unique sample of identical twins are compatible with the amyloid cascade hypothesis and thereby provide important new knowledge for clinical trial designs.

Plasma biomarkers predict amyloid pathology in cognitively normal monozygotic twins after 10 years.

Blood-based biomarkers could prove useful to predict Alzheimer's disease core pathologies in advance of clinical symptoms. Implementation of such biomarkers requires a solid understanding of their long-term dynamics and the contribution of confounding to their association with Alzheimer's disease pathology. Here we assess the value of plasma amyloid-ß1-42/1-40, phosphorylated-tau181 and glial fibrillary acidic protein to detect early Alzheimer's disease pathology, accounting for confounding by genetic and early environmental factors. Participants were 200 monozygotic twins, aged ≥60 years with normal cognition from the european medical information framework for Alzheimer's disease study. All twins had amyloid-ß status and plasma samples available at study enrolment. For 80 twins, additional plasma samples were available that had been collected approximately 10 years prior to amyloid-ß status assessment. Single-molecule array assays were applied to measure amyloid-ß1-42/1-40, phosphorylated-tau181 and glial fibrillary acidic protein. Predictive value of and longitudinal change in these biomarkers were assessed using receiver operating characteristic curve analysis and linear mixed models. Amyloid pathology could be predicted using blood-based biomarkers obtained at the time of amyloid status assessment (amyloid-ß1-42/1-40: area under the curve = 0.65, P = 0.01; phosphorylated-tau181: area under the curve = 0.84, P < 0.001; glial fibrillary acidic protein: area under the curve = 0.74, P < 0.001), as well as using those obtained 10 years prior to amyloid status assessment (amyloid-ß1-42/1-40: area under the curve = 0.69, P = 0.03; phosphorylated-tau181: area under the curve = 0.92, P < 0.001; glial fibrillary acidic protein: area under the curve = 0.84, P < 0.001). Longitudinally, amyloid-ß1-42/1-40 levels decreased [ß (SE) = -0.12 (0.01), P < 0.001] and phosphorylated-tau181 levels increased [ß (SE) = 0.02 (0.01), P = 0.004]. Amyloid-ß-positive individuals showed a steeper increase in phosphorylated-tau181 compared with amyloid-ß-negative individuals [ß (SE) = 0.06 (0.02), P = 0.004]. Also amyloid-ß-positive individuals tended to show a steeper increase in glial fibrillary acidic protein [ß (SE) = 0.04 (0.02), P = 0.07]. Within monozygotic twin pairs, those with higher plasma phosphorylated-tau181 and lower amyloid-ß1-42/1-40 levels were more likely to be amyloid-ß positive [ß (SE) = 0.95 (0.26), P < 0.001; ß (SE) = -0.28 (0.14), P < 0.05] indicating minimal contribution of confounding by genetic and early environmental factors. Our data support the use of amyloid-ß1-42/1-40, phosphorylated-tau181 and glial fibrillary acidic protein as screening tools for Alzheimer's disease pathology in the normal aging population, which is of importance for enrolment of high-risk subjects in secondary, or even primary, prevention trials. Furthermore, these markers show potential as low-invasive monitoring tool of disease progression and possibly treatment effects in clinical trials.

DNA methylation analysis is used to identify novel genetic loci associated with circulating fibrinogen levels in blood.

BACKGROUND: Fibrinogen plays an essential role in blood coagulation and inflammation. Circulating fibrinogen levels may be determined based on interindividual differences in DNA methylation at cytosine-phosphate-guanine (CpG) sites and vice versa. OBJECTIVES: To perform an EWAS to examine an association between blood DNA methylation levels and circulating fibrinogen levels to better understand its biological and pathophysiological actions. METHODS: We performed an epigenome-wide association study of circulating fibrinogen levels in 18 037 White, Black, American Indian, and Hispanic participants, representing 14 studies from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium. Circulating leukocyte DNA methylation was measured using the Illumina 450K array in 12 904 participants and using the EPIC array in 5133 participants. In each study, an epigenome-wide association study of fibrinogen was performed using linear mixed models adjusted for potential confounders. Study-specific results were combined using array-specific meta-analysis, followed by cross-replication of epigenome-wide significant associations. We compared models with and without CRP adjustment to examine the role of inflammation. RESULTS: We identified 208 and 87 significant CpG sites associated with fibrinogen levels from the 450K (p < 1.03 × 10-7) and EPIC arrays (p < 5.78 × 10-8), respectively. There were 78 associations from the 450K array that replicated in the EPIC array and 26 vice versa. After accounting for overlapping sites, there were 83 replicated CpG sites located in 61 loci, of which only 4 have been previously reported for fibrinogen. The examples of genes located near these CpG sites were SOCS3 and AIM2, which are involved in inflammatory pathways. The associations of all 83 replicated CpG sites were attenuated after CRP adjustment, although many remained significant. CONCLUSION: We identified 83 CpG sites associated with circulating fibrinogen levels. These associations are partially driven by inflammatory pathways shared by both fibrinogen and CRP.

Cardiorespiratory fitness, regular physical activity, and autonomic nervous system reactivity to laboratory and daily life stress.

The cross-stressor adaptation hypothesis-which posits that adjustment to physical stress as a result of regular physical activity and its effects on fitness crosses over to psychological stress reactivity-has been around for over four decades. However, the literature has been plagued by heterogeneities preventing definitive conclusions. We address these heterogeneity issues in a combined laboratory and daily life study of 116 young adults (M = 22.48 SD = 3.56, 57.76% female). The exposure, i.e., the potential driver of adaptation, was defined in three ways. First, a submaximal test was performed to obtain aerobic fitness measured as the VO2 max (kg/ml/min). Second, leisure time exercise behavior, and third, overall moderate-to-vigorous physical activity (MVPA), were obtained from a structured interview. Outcomes were autonomic nervous system (ANS) reactivity and affective responsiveness to stressors. ANS activity was measured continuously and expressed as inter-beat-interval (IBI), pre-ejection-period (PEP), respiratory sinus arrythmia (RSA), and non-specific Skin Conductance Responses (ns.SCR). Negative and positive affect were recorded after each experimental condition in the laboratory and hourly in daily life with a nine-item digital questionnaire. Linear regressions were performed between the three exposure measures as predictors and the various laboratory and daily life stress measurements as outcomes. Our results support the resting heart rate reducing effect of aerobic fitness and total MVPA in both the laboratory and daily life. We did not find evidence for the cross-stressor adaptation hypothesis, irrespective of ANS or affective outcome measure or whether the exposure was defined as exercise/MVPA or aerobic fitness.

Sociodemographic, lifestyle and clinical characteristics of energy-related depression symptoms: A pooled analysis of 13,965 depressed cases in 8 Dutch cohorts.

BACKGROUND: In a substantial subgroup of depressed patients, atypical, energy-related depression symptoms (e.g. increased appetite/weight, hypersomnia, loss of energy) tend to cluster with immuno-metabolic dysregulations (e.g. increased BMI and inflammatory markers). This clustering is proposed to reflect a more homogeneous depression pathology. This study examines to what extent energy-related symptoms are associated and share sociodemographic, lifestyle and clinical characteristics. METHODS: Data were available from 13,965 participants from eight Dutch cohorts with DSM-5 lifetime major depression assessed by the Lifetime Depression Assessment Self-report (LIDAS) questionnaire. Information on four energy-related depression symptoms were extracted: energy loss, increased appetite, increased weight, and hypersomnia. Tetrachoric correlations between these symptoms, and associations of these symptoms with sociodemographic (sex, age, education), lifestyle (physical activity, BMI, smoking) and clinical characteristics (age of onset, episode duration, history, treatment and recency, and self-reported comorbidity) were computed. RESULTS: Correlations between energy-related symptoms were overall higher than those with other depression symptoms and varied from 0.90 (increased appetite vs increased weight) to 0.11 (increased appetite vs energy loss). All energy-related symptoms were strongly associated with higher BMI and a more severe clinical profile. Patients with increased appetite were more often smokers, and only patients with increased appetite or weight more often had a self-reported diagnosis of PTSD (OR = 1.17, p = 2.91E-08) and eating disorder (OR = 1.40, p = 4.08E-17). CONCLUSIONS: The symptom-specific associations may have consequences for a profile integrating these symptoms, which can be used to reflect immuno-metabolic depression. They indicate the need to study immuno-metabolic depression at individual symptom resolution as a starting point.

Genetic Pathways Underlying Individual Differences in Regular Physical Activity.

Twin and family studies show a strong contribution of genetic factors to physical activity (PA) assessed by either self-report or accelerometers. PA heritability is around 43% across the lifespan. Genome-wide association studies have implied biological pathways related to exercise ability and enjoyment. A polygenic score based on genetic variants influencing PA could help improve the success of intervention programs.

Amyloid-ß and APOE genotype predict memory decline in cognitively unimpaired older individuals independently of Alzheimer's disease polygenic risk score.

BACKGROUND: What combination of risk factors for Alzheimer's disease (AD) are most predictive of cognitive decline in cognitively unimpaired individuals remains largely unclear. We studied associations between APOE genotype, AD-Polygenic Risk Scores (AD-PRS), amyloid-ß pathology and decline in cognitive functioning over time in a large sample of cognitively unimpaired older individuals. METHODS: We included 276 cognitively unimpaired older individuals (75 ± 10 years, 63% female) from the EMIF-AD PreclinAD cohort. An AD-PRS was calculated including 83 genome-wide significant variants. The APOE gene was not included in the PRS and was analyzed separately. Baseline amyloid-ß status was assessed by visual read of [18F]flutemetamol-PET standardized uptake value images. At baseline and follow-up (2.0 ± 0.4 years), the cognitive domains of memory, attention, executive function, and language were measured. We used generalized estimating equations corrected for age, sex and center to examine associations between APOE genotype and AD-PRS with amyloid-ß status. Linear mixed models corrected for age, sex, center and education were used to examine associations between APOE genotype, AD-PRS and amyloid-ß status, and their interaction on changes in cognitive functioning over time. RESULTS: Fifty-two participants (19%) had abnormal amyloid-ß, and 84 participants (31%) carried at least one APOE ε4 allele. APOE genotype and AD-PRS were both associated with abnormal amyloid-ß status. Increasingly more risk-full APOE genotype, a high AD-PRS and an abnormal amyloid-ß status were associated with steeper decline in memory functioning in separate models (all p ≤ 0.02). A model including 4-way interaction term (APOE×AD-PRS×amyloid-ß×time) was not significant. When modelled together, both APOE genotype and AD-PRS predicted steeper decline in memory functioning (APOE ß(SE)=-0.05(0.02); AD-PRS ß(SE)=-0.04(0.01)). Additionally, when modelled together, both amyloid-ß status and AD-PRS predicted a steeper decline in memory functioning (amyloid-ß ß(SE)=-0.07(0.04); AD-PRS ß(SE)=-0.04(0.01)). Modelling both APOE genotype and amyloid-ß status, we observed an interaction, in which APOE genotype was related to steeper decline in memory and language functioning in amyloid-ß abnormal individuals only (ß(SE)=-0.13(0.06); ß(SE)=-0.22(0.07), respectively). CONCLUSION: Our results suggest that APOE genotype is related to steeper decline in memory and language functioning in individuals with abnormal amyloid-ß only. Furthermore, independent of amyloid-ß status other genetic risk variants contribute to memory decline in initially cognitively unimpaired older individuals.